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We characterize the synchronization of an array of coupled chaotic elements as a phase transition where order 
parameters related to the joint probability at two sites obey power laws versus the mutual coupling strength; the 
phase transition corresponds to a change in the exponent of the power law. Since these studies are motivated 
by the behaviour of the cortical neurons in cognitive tasks, we account for the short time available to any brain 
decision by studying how the mutual coupling affects the transient behaviour of a synchronization transition 
over a fixed time interval. We present a novel feature, namely, the absence of decay of the initial defect density 
for small coupling. 
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Temporal coding vs. rate coding for the neural-based infor- 
mation has been open to debate in the neuroscience literature 
1 1 ] . In rate coding only the mean frequency of spikes over a 
time interval matters, thus requiring a suitable counting inter- 
val, which seems unfit for fast decision tasks. Temporal cod- 
ing assigns importance to the precise timing and coordination 
of action potentials or "spikes". A special type of temporal 
coding is synchrony, whereby information is encoded by the 
synchronous firing of action potentials of all neurons of a cor- 
tical module. 

Let a laboratory animal be exposed to a visual field con- 
taining two separate objects, both made of the same visual 
elements. Since each receptive field isolates a specific detail 
of an object, we should expect a corresponding large set of 
individual responses. On the contrary, all the cortical neu- 
rons whose receptive fields are pointing to a specific object 
synchronize their spikes, and as a consequence the visual cor- 
tex organizes into separate neuron groups oscillating on two 
distinct spike trains for the two objects (feature binding). Ex- 
perimental evidence was obtained by insertion in the cortical 
tissue of animals of microelectrodes each sensing a single neu- 
ron (2l • Indirect evidence of synchronization has been reached 
for human subjects as well, by processing the EEG (electro- 
encephalo-gram) data 

The dynamics ruling the above facts should be described 
by a convenient dynamical model, independently of the bio- 
logical components. Since feature binding results from the 
readjustment of the temporal positions of the spikes, the most 
convenient strategy for this is to rely on the properties of suit- 
able chaotic oscillators (c-o) 01 • 

In this paper we explore the conditions under which an as- 
sembly of coupled c-o reaches mutual synchronization, thus 
displaying a coherent behaviour. We compare two dynamical 
models, namely HC |5] and Roessler |6]. Homoclinic chaos 
(HC) |5|] appears as the optimal strategy for a time code shared 
by a large crowd of identical coupled objects. Indeed HC pro- 
vides at each pseudo-period (or inter spike interval = ISI) 
the alternation of a regular large spike (a) and a small chaotic 
background (b). (b) is the sensitive region where the activation 



from neighbors occurs, while the spike (a) provides a suitable 
signal to activate the coupling. Whence, a chain of weakly 
coupled c-o of this kind will easily synchronize, reaching a 
state common to all sites. 

We consider an array of identical c-o with nearest neigh- 
bour bidirectional coupling, ruled by L2J 
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The index i denotes the ith site po- 



. , N), and dots denote temporal derivatives. 
The model |7] was first introduced to describe the chaotic dy- 
namics of a molecular (CO2) laser with feedback. For each 
site i, variable xi represents the laser intensity, X2 the popu- 
lation inversion between the two levels resonant with the ra- 
diation field, and xq the feedback voltage which controls the 
cavity losses. The auxiliary variables x^, X4 and X5 account 
for molecular exchanges within the CO2 molecule. We con- 
sider identical systems at each site i; as for the parameters, 
their physical meaning has been discussed elsewhere. Their 
values are: fco = 28.5714, ki = 4.5556, 71 = 10.0643, 
72 = 1.0643, g = 0.05, po = 0.016, z = 10, /3 = 0.4286, 
a = 32.8767, r = 160 and bo = 0.1032. The mutual cou- 
pling consists of adding to the equation for on each site 
a function of the intensity of the neighboring oscillators. 
The term {x\) represents the average value of the x\ variable, 
calculated as a moving average over a long time. The system 
is integrated by means of Bulirsch-Stoer predictor-corrector 
method with open boundary conditions. 

Any c-o where chaos is due to the homoclinic (or hetero- 
clinic) return to a saddle focus displays a high sensitivity to 
an external perturbation just in the neighborhood of the saddle 
|5]. We use the generic attribution of HC for systems which 
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have a saddle focus instability f^. In view of this high sen- 
sitivity, we expect that they synchronize not only under an 
external driving |9|, but also for a convenient mutual coupling 
strength. A quantitative indicator of this sensitivity is rep- 
resented by the so-called propensity to synchronization |1(|] 
whereby different chaotic systems can be compared. Further- 
more, for an array of coupled c-o, the onset of a collective 
synchronization was explored numerically for some selected 
values of the coupling strength e, showing that the lack of 
synchronization manifests itself as phase slips which appear 
as dislocations in a space-time plot (space being the axis of 
site positions, time being the point-like occurrence of a spike) 
fl ill . Spacewise, a tiny variation of e changes dramatically 
the size of the synchronized domain 1 10]. Here we study in 
a systematic way the onset of synchronization. Introducing 
a suitable order paramater, we prove that it has a power law 
dependence on e (control parameter); thus characterizing syn- 
chronization as a phase transition. 

The most prominent feature of HC is that the time signal 
consists of almost identical spikes occuring at chaotic times, 
with an {I SI) (average interspike interval) much larger than 
the spike duration. Synchronization of two adjacent sites im- 
plies that pairs of spikes occuring on the two sites be separated 
by much less then {I SI) . Precisely, when the spikes are sep- 
arated by less then 10%{ISI) we will say that the two sites 
are locally synchronized and attribute a mutual probability of 
spike occurance Pi.i+i = 1. If the separation is larger then 
10%{ISI) then we take Pi^i+i — and claim that locally 
synchronization is lost, because of a phase slip |9]. Based on 
this set of joint probabilities we build a mutual information: 
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Let us consider a linear array of N nearest neighbours coupled 
sites. As we increase the coupling parameter e, the /y vs the 
separation \j ~i\ — x goes as Fig.^ We plot the average gra- 
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FIG. 1 : Mutual information of two sites i and j versus the site sepa- 
ration X — |i — j| for various amounts e of the coupling strength. 
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of lij versus e (Fig.|2p). While initially goes as 
e~7,ate=:ei there is a sharp transition point to a new power 
law dependence as . Eventually at e = £2 we have a sec- 
ond transition to e^^. In order to give a meaning to the two 
transition points ei.2, we plot the "cluster size" vs e. Using 



the method for cluster counting described in Ref. fT?], two 
sites belong to the same cluster if 1^ > 9 where 6* is a fixed 
threshold. This way, we define "cluster size" as the size of the 
largest cluster which occurs within the array during the time 
evolution. The cluster size increases initially as , then as 
between ei and e2 and finally as above e = £2 (Fig. |2j)). 
Now, an increase as e 2 is typical of a linear diffusive process. 
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FIG. 2: Linear array of 40 coupled HC c-o. (a) Average gradient 
of the cluster size A*'. It coincides with the gradient of the mutual 
information (see Fig.0 averaged over all sites. The two cusp points 
are ei = 0.0245 and £2 ~ 0.049. (b) As in a). Average cluster size 
versus e. The two cusp points are ei = 0.0245 and e2 = 0.049 (the 
second one corresponding to saturation). (c) Array of 40 Roessler c-o. 
The cluster size gradient displays a single cusp at ei = 0.00405. (d) 
As in c). Cluster size versus e. There is only one cusp corresponding 
to saturation at ei = 0.00405. At variance with HC there is no 
intermediate speed up. 

Indeed, the correlation domain in linear diffusion increases as 
Ic oc \/ei, thus, for a fixed observation time t, we have the 
power law Ic (x . We therefore consider the steep increase 
beyond the transition point e > £2 as a nontrivial nonlinear 
effect. We take the change of slope at ei as the onset of syn- 
chronization. We wish to show that this steep synchronization 
is a peculiarity of HC systems. For this purpose, we study an 
array of chaotic oscillators where no spike, but rather phase 
synchronization occurs. Such is the case of coupled Roessler 
oscillators |6], where the chaotic oscillations are character- 
ized by a main oscillation frequency, and chaos appears as a 
spread in the amplitude of different cycles; here in fact the 
array tends toward a phase synchronization. The dynamics is 
ruled by O 
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where a — 0.15, b = 0.2 and c = 10. Independently from the 
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chaotic amplitude, the phase difference at two different sites 
can vary smoothly from to vr. Calling T ^ ^ the period of 
the main frequency ujq, then the joint probability of two sites 
can be taken as: 
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and it goes with continuity from 1 (Ty = 0) to (T^ — _._ 2 , 
Using such a probability we evaluate the mutual information 
and the cluster size. Fig. ^ and Fig. |2ll. They show that £2 



(a) 



(b) 



^ 10 
i 



— ^ ■-■ — ■ 




FIG. 3: Spatiotemporal evolutions of HC systems for coupling con- 
stant (a) e = 0.03 and (b) e — 0.1; time in the vertical axis is given 
in (ISI) units. 

is missing, and above ei the average gradient changes from 
to e^i. Notice that, up to total invasion of the available 
domain, the progress of phase synchronization grows as e^, 
as in linear diffusion. 

This fact confirms that the steep rise of synchronization be- 
yond 62 is a peculiar aspect of HC, or, more generally, of any 
chaotic dynamical system where chaos is due to the homo- 
clinic return to a saddle focus. We use the generic attribution 
of HC for a large class of systems, with a saddle focus insta- 
bility 1 8]. This includes a class B laser with feedback |7], the 
Hodgkin-Huxley model for the production of action potentials 
in a neuron membrane and the Hindmarsh-Rose model of 
generation of spike bursts [15]. As a partial conclusion, HC is 
a very plausible model to explain feature binding in an array 
of neurons belonging to the same cortical module. 

The onset of synchronization, either spike in HC or phase 
in Roessler, is described by order parameters (the gradient of 
the mutual information and the cluster size) which undergo a 
phase transition as functions of a control parameter (the mu- 
tual coupling strength). The signature of the phase transition 
is the change of the power law behaviour, evaluated for long 
time well beyond any initial transient. 

However, a perceptual task must be completed within a 
short time after the arrival of external stimuli, in order for the 
living organism to take vital decisions as a response to those 
stimuli. This time for human subjects is around 200 ms which 
corresponds to about 10(75'/) in the so called 7 band which 
provides relevant cues in visual tasks |6]. 

We then consider a transient phase transition, requiring that 
all sites be synchronized within less than 20 {I SI). Notice that 
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FIG. 4: Transient times versus cluster size for (a) HC and (b) 
Roessler c-o calculated for selected values of coupling constant dur- 
ing the limited time period taken as approximately 20 spikes. Tran- 
sient time is represented in the units of the average interspike interval 
{ISI) in the case of HC and in units of constant period T in the case 
of Roessler. 
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FIG. 5: Linear array of 40 coupled HC (a) and Roessler (b) c-o. 
Average gradient of the transient time ttrans versus coupling e. The 
two cusp points in (a) are at ei — 0.0245 and €2 — 0.049 and the 
cusp in (b) is at ei = 0.00405. The cusp points in both cases coincide 
with those of Fig.|2| 



the criterion given above for = 1 in HC does not iden- 

tify synchronization with isochronism, but it just rules out a 
phase slip (one spike more or less) between site i and j. In a 
space-time plot where spike occurrence is denoted by a dot, 
a phase slip would appear as a dislocation (lack of connect- 
edness) in an a rray of lines joining the spike occurrences at 
different sites llOll . Space-time plots are shown in Fig. [3^ 
(e = 0.03) and Fig.lSj) (e = 0.1) where the time is given in 
(ISI) units; the plots show time intervals starting from t — 
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as well as from t = 500. 

For e = 0.03, the average number of defects does not decay 
even for long times; indeed a plot of the same interval starting 
from t — 500 displays the same number of defects on average; 
defects die out at one site and are born again at another one in 
course of time, due to the temporal relations between adjacent 
sites necessary to induce the escape from the saddle region 
(see the detailed discussion for unidirectional coupling in Ref. 
llal). On the contrary, for e = 0.1, defects decays within 
20(/S'/). 

In Fig. 13 we show the transient time (in {I SI) units) nec- 
essary to reach a cluster size x (x can vary from 1 to 40), 
calculated for an ensemble of random initial conditions. The 
apparent monotonic increase with time of the cluster size for 
low e is an ensemble averaging contrivance, irrelevant for the 
single trial. Having the distribution of transient times over 
the cluster sizes for different couplings we consider the aver- 
aged space gradient ^MhzEzsl^ This quantity yields informa- 
tion on how the transient times needed to reach a given cluster 
size changes with the coupling and provides a global picture 
of the transient dynamics in the spatiotemporal system. For 
low e we have a flat dependence on e which means that the 
coupling is not effective at all in inducing synchronization, 
neither in HC (Fig.|5t) nor in Roessler (Fig.lSji); the critical 
values ei = 0.0245 and 62 = 0.049 sign the starting point 
for a power law decay of that gradient in HC; similarly does 
ei = 0.00405 for Roessler 

Let us elaborate on the persistence of initial defects and 
show that it is peculiar of synchronized HC arrays. In HC a de- 
fect corresponds to a sudden jump of pij from 1 to or vicev- 
ersa. In order to investigate this new feature in another system, 
we introduce a symbolic dynamics in Roessler as well. In Eq. 
Elpij goes smoothly from 1 at Tij = to at = ^; we 
replace that smooth behaviour with a jumpy one, by taking 
Pij = 1 when Tij > i and pij — when Tij < i. This 
way, we obtain space-time plots of defects for Roessler which 
resemble those of Fig.|3] However, here defects always decay 
for long times, as checked for e = 0.001 and e = 0.0001, 
and thus the array is fully synchronized after a finite time, 
even for very small e. The flat region below ei = 0.00405 in 
Fig-Eb is due to the fact that we were studying the transition 
to synchronization for short times, up to 20r. Thus, there ex- 
ists a crucial difference between HC and Roessler: for small 
e in HC, even after reaching the synchronization of the whole 
array, the defects may reappear and then again disappear, on 
the contrary in Roessler, once full synchronization is reached, 
the defects do not reappear Transient defects considered in 
sudden symmetry changes of Ginzburg-Landau type systems 
decay always, as shown theoretically ll^ and verified exper- 
imentally in He3 1 18j and in nonlinear optics fl9ll . Instead, 
below synchronization threshold of an HC array, there is a 
steady defect sea. It would be interesting to explore the ap- 
pearance of such a feature in situations different from that here 
considered. Furthermore, this report has been limited to spon- 
taneous synchronization; from a biological point of view it is 



crucial to explore the general scaling laws in presence of ex- 
ternal forcing as already done, but only for some numerical 
examples II ill . 
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parameters related to the joint probability at two sites obey power laws versus the mutual coupling strength; the 
phase transition corresponds to a change in the exponent of the power law. Since these studies are motivated 
by the behaviour of the cortical neurons in cognitive tasks, we account for the short time available to any brain 
decision by studying how the mutual coupling affects the transient behaviour of a synchronization transition 
over a fixed time interval. We present a novel feature, namely, the absence of decay of the initial defect density 
for small coupling. 
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Temporal coding vs. rate coding for the neural-based infor- 
mation has been open to debate in the neuroscience literature 
1 1 ] . In rate coding only the mean frequency of spikes over a 
time interval matters, thus requiring a suitable counting inter- 
val, which seems unfit for fast decision tasks. Temporal cod- 
ing assigns importance to the precise timing and coordination 
of action potentials or "spikes". A special type of temporal 
coding is synchrony, whereby information is encoded by the 
synchronous firing of action potentials of all neurons of a cor- 
tical module. 

Let a laboratory animal be exposed to a visual field con- 
taining two separate objects, both made of the same visual 
elements. Since each receptive field isolates a specific detail 
of an object, we should expect a corresponding large set of 
individual responses. On the contrary, all the cortical neu- 
rons whose receptive fields are pointing to a specific object 
synchronize their spikes, and as a consequence the visual cor- 
tex organizes into separate neuron groups oscillating on two 
distinct spike trains for the two objects (feature binding). Ex- 
perimental evidence was obtained by insertion in the corti- 
cal tissue of animals of microelectrodes each sensing a single 
neuron j^|3|]. Indirect evidence of synchronization has been 
reached for human subjects as well, by processing the EEG 
(electro-encephalo-gram) data 

The dynamics ruling the above facts should be described 
by a convenient dynamical model, independently of the bio- 
logical components. Since feature binding results from the 
readjustment of the temporal positions of the spikes, the most 
convenient strategy for this is to rely on the properties of suit- 
able chaotic oscillators (c-o) 01 • 

In this paper we explore the conditions under which an as- 
sembly of coupled c-o reaches mutual synchronization, thus 
displaying a coherent behaviour. We compare two dynamical 
models, namely HC |5] and Roessler |6]. Homoclinic chaos 
(HC) |5|] appears as the optimal strategy for a time code shared 
by a large crowd of identical coupled objects. Indeed HC pro- 
vides at each pseudo-period (or inter spike interval = ISI) 
the alternation of a regular large spike (a) and a small chaotic 
background (b). (b) is the sensitive region where the activation 



from neighbors occurs, while the spike (a) provides a suitable 
signal to activate the coupling. Whence, a chain of weakly 
coupled c-o of this kind will easily synchronize, reaching a 
state common to all sites. 

We consider an array of identical c-o 0] with nearest neigh- 
bour bidirectional coupling, ruled by 
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The index i denotes the ith site po- 



. , N), and dots denote temporal derivatives. 
The model |7] was first introduced to describe the chaotic dy- 
namics of a molecular (CO2) laser with feedback. For each 
site i, variable xi represents the laser intensity, X2 the popu- 
lation inversion between the two levels resonant with the ra- 
diation field, and xq the feedback voltage which controls the 
cavity losses. The auxiliary variables x^, X4 and X5 account 
for molecular exchanges within the CO2 molecule. We con- 
sider identical systems at each site i; as for the parameters, 
their physical meaning has been discussed elsewhere. Their 
values are: fco = 28.5714, ki = 4.5556, 71 = 10.0643, 
72 = 1.0643, g = 0.05, po = 0.016, z = 10, /3 = 0.4286, 
a = 32.8767, r = 160 and bo = 0.1032. The mutual cou- 
pling consists of adding to the equation for on each site 
a function of the intensity of the neighboring oscillators. 
The term {x\) represents the average value of the x\ variable, 
calculated as a moving average over a long time. The system 
is integrated by means of Bulirsch-Stoer predictor-corrector 
method with open boundary conditions. 

Any c-o where chaos is due to the homoclinic (or hetero- 
clinic) return to a saddle focus displays a high sensitivity to 
an external perturbation just in the neighborhood of the saddle 
|5]. We use the generic attribution of HC for systems which 
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have a saddle focus instability f^. In view of this high sen- 
sitivity, we expect that they synchronize not only under an 
external driving |9|, but also for a convenient mutual coupling 
strength. A quantitative indicator of this sensitivity is rep- 
resented by the so-called propensity to synchronization |1(|] 
whereby different chaotic systems can be compared. Further- 
more, for an array of coupled c-o, the onset of a collective 
synchronization was explored numerically for some selected 
values of the coupling strength e, showing that the lack of 
synchronization manifests itself as phase slips which appear 
as dislocations in a space-time plot (space being the axis of 
site positions, time being the point-like occurrence of a spike) 
fl ill . Spacewise, a tiny variation of e changes dramatically 
the size of the synchronized domain 1 10]. Here we study in 
a systematic way the onset of synchronization. Introducing 
a suitable order paramater, we prove that it has a power law 
dependence on e (control parameter); thus characterizing syn- 
chronization as a phase transition. 

The most prominent feature of HC is that the time signal 
consists of almost identical spikes occuring at chaotic times, 
with an {I SI) (average interspike interval) much larger than 
the spike duration. Synchronization of two adjacent sites im- 
plies that pairs of spikes occuring on the two sites be separated 
by much less then {I SI) . Precisely, when the spikes are sep- 
arated by less then 10%{ISI) we will say that the two sites 
are locally synchronized and attribute a mutual probability of 
spike occurance Pi.i+i = 1. If the separation is larger then 
10%{ISI) then we take Pi^i+i — and claim that locally 
synchronization is lost, because of a phase slip |9]. Based on 
this set of joint probabilities we build a mutual information: 
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Let us consider a linear array of N nearest neighbours coupled 
sites. As we increase the coupling parameter e, the /y vs the 
separation \j ~i\ — x goes as Fig.^ We plot the average gra- 
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FIG. 1 : Mutual information of two sites i and j versus the site sepa- 
ration X — |i — j| for various amounts e of the coupling strength. 
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of lij versus e (Fig.|2p). While initially goes as 
e~7,ate=:ei there is a sharp transition point to a new power 
law dependence as . Eventually at e = £2 we have a sec- 
ond transition to e^^. In order to give a meaning to the two 
transition points ei.2, we plot the "cluster size" vs e. Using 



the method for cluster counting described in Ref. fT?], two 
sites belong to the same cluster if 1^ > 9 where 6* is a fixed 
threshold. This way, we define "cluster size" as the size of the 
largest cluster which occurs within the array during the time 
evolution. The cluster size increases initially as , then as 
between ei and e2 and finally as above e = £2 (Fig. |2j)). 
Now, an increase as e 2 is typical of a linear diffusive process. 
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FIG. 2: Linear array of 40 coupled HC c-o. (a) Average gradient 
of the cluster size A*'. It coincides with the gradient of the mutual 
information (see Fig.0 averaged over all sites. The two cusp points 
are ei = 0.0245 and £2 ~ 0.049. (b) As in a). Average cluster size 
versus e. The two cusp points are ei = 0.0245 and e2 = 0.049 (the 
second one corresponding to saturation). (c) Array of 40 Roessler c-o. 
The cluster size gradient displays a single cusp at ei = 0.00405. (d) 
As in c). Cluster size versus e. There is only one cusp corresponding 
to saturation at ei = 0.00405. At variance with HC there is no 
intermediate speed up. 

Indeed, the correlation domain in linear diffusion increases as 
Ic oc \/ei, thus, for a fixed observation time t, we have the 
power law Ic (x . We therefore consider the steep increase 
beyond the transition point e > £2 as a nontrivial nonlinear 
effect. We take the change of slope at ei as the onset of syn- 
chronization. We wish to show that this steep synchronization 
is a peculiarity of HC systems. For this purpose, we study an 
array of chaotic oscillators where no spike, but rather phase 
synchronization occurs. Such is the case of coupled Roessler 
oscillators |6], where the chaotic oscillations are character- 
ized by a main oscillation frequency, and chaos appears as a 
spread in the amplitude of different cycles; here in fact the 
array tends toward a phase synchronization. The dynamics is 
ruled by O 
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±1 = b + xl{x\-c) 



(3) 



where a — 0.15, b = 0.2 and c = 10. Independently from the 
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chaotic amplitude, the phase difference at two different sites 
can vary smoothly from to vr. Calling T ^ ^ the period of 
the main frequency ujq, then the joint probability of two sites 
can be taken as: 



(a) 



(b) 



T 



(4) 



and it goes with continuity from 1 (Ty = 0) to (T^ — _._ 2 , 
Using such a probability we evaluate the mutual information 
and the cluster size. Fig. ^ and Fig. |2ll. They show that £2 
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FIG. 3: Spatiotemporal evolutions of HC systems for coupling con- 
stant (a) e = 0.03 and (b) e — 0.1; time in the vertical axis is given 
in (ISI) units. 

is missing, and above ei the average gradient changes from 
to e^i. Notice that, up to total invasion of the available 
domain, the progress of phase synchronization grows as e^, 
as in linear diffusion. 

This fact confirms that the steep rise of synchronization be- 
yond 62 is a peculiar aspect of HC, or, more generally, of any 
chaotic dynamical system where chaos is due to the homo- 
clinic return to a saddle focus. We use the generic attribution 
of HC for a large class of systems, with a saddle focus insta- 
bility 1 8]. This includes a class B laser with feedback |7], the 
Hodgkin-Huxley model for the production of action potentials 
in a neuron membrane and the Hindmarsh-Rose model of 
generation of spike bursts [15]. As a partial conclusion, HC is 
a very plausible model to explain feature binding in an array 
of neurons belonging to the same cortical module. 

The onset of synchronization, either spike in HC or phase 
in Roessler, is described by order parameters (the gradient of 
the mutual information and the cluster size) which undergo a 
phase transition as functions of a control parameter (the mu- 
tual coupling strength). The signature of the phase transition 
is the change of the power law behaviour, evaluated for long 
time well beyond any initial transient. 

However, a perceptual task must be completed within a 
short time after the arrival of external stimuli, in order for the 
living organism to take vital decisions as a response to those 
stimuli. This time for human subjects is around 200 ms which 
corresponds to about 10(75'/) in the so called 7 band which 
provides relevant cues in visual tasks |6]. 

We then consider a transient phase transition, requiring that 
all sites be synchronized within less than 20 {I SI). Notice that 
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FIG. 4: Transient times versus cluster size for (a) HC and (b) 
Roessler c-o calculated for selected values of coupling constant dur- 
ing the limited time period taken as approximately 20 spikes. Tran- 
sient time is represented in the units of the average interspike interval 
{ISI) in the case of HC and in units of constant period T in the case 
of Roessler. 
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FIG. 5: Linear array of 40 coupled HC (a) and Roessler (b) c-o. 
Average gradient of the transient time ttrans versus coupling e. The 
two cusp points in (a) are at ei — 0.0245 and €2 — 0.049 and the 
cusp in (b) is at ei = 0.00405. The cusp points in both cases coincide 
with those of Fig.|2| 



the criterion given above for = 1 in HC does not iden- 

tify synchronization with isochronism, but it just rules out a 
phase slip (one spike more or less) between site i and j. In a 
space-time plot where spike occurrence is denoted by a dot, 
a phase slip would appear as a dislocation (lack of connect- 
edness) in an a rray of lines joining the spike occurrences at 
different sites llOll . Space-time plots are shown in Fig. [3^ 
(e = 0.03) and Fig.lSj) (e = 0.1) where the time is given in 
(ISI) units; the plots show time intervals starting from t — 



4 



as well as from t = 500. 

For e = 0.03, the average number of defects does not decay 
even for long times; indeed a plot of the same interval starting 
from t — 500 displays the same number of defects on average; 
defects die out at one site and are born again at another one in 
course of time, due to the temporal relations between adjacent 
sites necessary to induce the escape from the saddle region 
(see the detailed discussion for unidirectional coupling in Ref. 
llal). On the contrary, for e = 0.1, defects decays within 
20(/S'/). 

In Fig. 13 we show the transient time (in {I SI) units) nec- 
essary to reach a cluster size x (x can vary from 1 to 40), 
calculated for an ensemble of random initial conditions. The 
apparent monotonic increase with time of the cluster size for 
low e is an ensemble averaging contrivance, irrelevant for the 
single trial. Having the distribution of transient times over 
the cluster sizes for different couplings we consider the aver- 
aged space gradient ^MhzEzsl^ This quantity yields informa- 
tion on how the transient times needed to reach a given cluster 
size changes with the coupling and provides a global picture 
of the transient dynamics in the spatiotemporal system. For 
low e we have a flat dependence on e which means that the 
coupling is not effective at all in inducing synchronization, 
neither in HC (Fig.|5t) nor in Roessler (Fig.lSji); the critical 
values ei = 0.0245 and 62 = 0.049 sign the starting point 
for a power law decay of that gradient in HC; similarly does 
ei = 0.00405 for Roessler 

Let us elaborate on the persistence of initial defects and 
show that it is peculiar of synchronized HC arrays. In HC a de- 
fect corresponds to a sudden jump of pij from 1 to or vicev- 
ersa. In order to investigate this new feature in another system, 
we introduce a symbolic dynamics in Roessler as well. In Eq. 
Elpij goes smoothly from 1 at Tij = to at = ^; we 
replace that smooth behaviour with a jumpy one, by taking 
Pij = 1 when Tij > i and pij — when Tij < i. This 
way, we obtain space-time plots of defects for Roessler which 
resemble those of Fig.|3] However, here defects always decay 
for long times, as checked for e = 0.001 and e = 0.0001, 
and thus the array is fully synchronized after a finite time, 
even for very small e. The flat region below ei = 0.00405 in 
Fig-Eb is due to the fact that we were studying the transition 
to synchronization for short times, up to 20r. Thus, there ex- 
ists a crucial difference between HC and Roessler: for small 
e in HC, even after reaching the synchronization of the whole 
array, the defects may reappear and then again disappear, on 
the contrary in Roessler, once full synchronization is reached, 
the defects do not reappear Transient defects considered in 
sudden symmetry changes of Ginzburg-Landau type systems 
decay always, as shown theoretically ll^ and verified exper- 
imentally in He3 1 18j and in nonlinear optics fl9ll . Instead, 
below synchronization threshold of an HC array, there is a 
steady defect sea. It would be interesting to explore the ap- 
pearance of such a feature in situations different from that here 
considered. Furthermore, this report has been limited to spon- 
taneous synchronization; from a biological point of view it is 



crucial to explore the general scaling laws in presence of ex- 
ternal forcing as already done, but only for some numerical 
examples II ill . 
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